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We report the current status of the on-going lattice-QCD calculations of nucleon isovector ax- 
ial charge, gA, using the RBC/UKQCD 2+1-flavor dynamical domain-wall fermion ensembles 
at lattice cutoff of about a -1 = 1.4 GeV in a spatial volume (L = 4.6fm) 3 . The result from 
the ensemble with 11% = 250 MeV pion mass, corresponding to the finite-size scaling parameter 
m K L ~ 5.8, agrees well with an earlier result at a = 1.7 GeV, L = 2.8 fm, and m n = 420 MeV, 
with similar m K L. This suggests the systematic error from excited-state contamination is small 
in both ensembles and about 10-% deficit in g& we are observing is likely a finite-size effect that 
scales with m K L. We also report the result from the lighter, m K = 170 MeV ensemble. 
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1. Introduction 

RBC and UKQCD collaborations have been jointly calculating nucleon-structure observables 
such as form factors and low moments of structure functions [1, 2, 3, 4, 6, 5] using their 2+1- 
flavor dynamical Domain-Wall Fermions (DWF) lattice-QCD ensembles [7, 8, 9]. Here we report 
the ongoing calculation using the latest ensembles with pion mass, m n , as low as 250 and 170 
MeV with a linear spatial extent, L, of about 4.6 fm [9], corresponding to the finite-size scaling 
parameter, m % L, of about 6 and 4 respectively. 

Using our earlier ensembles with pion mass down to about 330 MeV in a linear spatial extent of 
about 2.8 fm, we observed a significant deficit in the isovector axial charge, gA, of nucleon, of about 
10% compared with the experiment [1]. By comparing with many other similar lattice numerical 
calculations with relatively heavier mass and smaller spatial extents that showed similar trends, 
we conjectured the deficit is likely a finite-size effect from relatively small values of the scaling 
parameter, m K L [1]. This can be interpreted as the first ever concrete evidence for the "virtual pion 
cloud" surrounding the nucleon: Such a cloud carries isovector axialvector current of the nucleon, 
and when the pion is light, can fluctuate far away from the nucleon center of mass resulting in a 
deficit if the volume is insufficient. However some questions against this interpretation were raised 
in the past couple of years, that the deficit we and many others had seen might have been systematic 
errors from excited- state contamination [10, 11]. 

The preliminary results of our ongoing calculations, with different gauge action, different 
lattice cut off, much lighter pion mass and much larger lattice spatial extent are already sufficient 
to exclude such excited-state contamination from at least two of our calculations: the earlier one at 
rriji = 420 MeV and L = 2.8 fm and the ongoing one at m n = 250 MeV and L = 4.6 fm that share 
about the same value of m n L = 5.8, and reinforce our conjecture of scaling with m % L. 

Details of form factors are reported separately by Meifeng Lin in these proceedings [12]. 

2. Numerics 

The joint RBC and UKQCD 2+1 -flavor dynamical DWF ensembles [9] we are using in this 
report are generated with a new gauge action, a combination of conventional Iwasaki rectangular- 
improved action [13] at the inverse-squared gauge coupling of /3 = 1.75 and a multiplicative 
disolocation-suppressing determinant ratio (DSDR) factor [14, 15, 16]. This combination allows us 
calculations with relatively low lattice cut-off momentum of aT x ~ 1.371(8) GeV without driving 
up the residual chiral symmetry breaking so the DWF residual mass is m res a ~ 0.002 while main- 
taining reasonable distribution in gauge field homotopy. Consequently the linear spatial extent of 
the lattice, with 32 lattice sites in one direction, is about 4.6 fm, and allows us to calculate at much 
closer to physical pion mass, m% ~ 170 and 250 MeV, or input quark mass, m uc ia = 0.001 and 
0.0042. ALCF, a BG/P facility, was used to generate the gauge configurations used in this report. 

We use standard ratios, C^(f s i n k,?)/C2 P t(4ink), of the nucleon three-point, C 3 ' t (?sink>0 = 
Ha,pr a p{N p (t sink )O(t)N a (0)), to two-point, C 2pt (f s ink) = La,/, (1/2) (1 + 7i) a p{N p (tank)N a (0)), 
functions to calculate the nucleon observables, O, [17]. We use N = e a b c {u^CYsdb)u c as our nu- 
cleon operator, Gaussian smearing [18, 19] for the source and point sink, and appropriate projec- 
tions r. As was reported earlier, we optimized our gauge-invariant Gaussian-smearing width and 
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Figure 1: Left: nucleon mass effective -mass plot for RBC/UKQCD (2+l)-fiavor dynamical ensembles with 
lattice cut off of a' 1 = 1.371(8) GeV and pion mass m n — 250 MeV (red) and 170 MeV (black). We estimate 
the nucleon mass as m N = 0.769(5) or 1.054(7) GeV for the former and m N = 0.718(6) or 0.984(8) GeV for 
the latter. Right: mjy plotted against nv\ from recent RBC+UKQCD 2+1-flavor dynamical DWF ensembles. 

source-sink separation, t s - m ^, so as to remove the excited-state contamination below our statistical 
error: we chose the Gaussian width of six, and sour-sink separation of nine [4, 5, 6]. We use 165 
configurations for the heavy, m n = 250 MeV, ensemble, every eighth trajectory from the trajectory 
608 to 1920, and 1 14 configurations for the light, 170-MeV, ensemble, from 508 to 1412. For this 
report we have four source positions at t = 0, 16, 32 and 48 for each configuration. The calculations 
are done on the RIKEN RICC computer cluster for the former, and NSF Teragrid/XSEDE clusters 
for the latter. We obtain nucleon mass estimates of m^/ = 0.769(5) or 1.054(7) GeV for the former 
and mjy = 0.718(6) or 0.984(8) GeV for the latter (see Fig. 1.) These values fall on a smooth line 
with our earlier results to the physical point. 

3. Vector and axial charges and their ratio 

We calculate the isovector vector charge, gv, of nucleon from the time component of the local 
current. It deviates from unity and gives the inverse of the non-perturbative renormalization, Zy, 
of the current: gv = 1/Zy, up to small 0(a 2 ) correction. As is shown in the left pane of Fig. 
2, we obtain very clean and accurate determination of gv = 1.450(4) for the heavy, 250-MeV, 
ensemble and 1.447(9) for the light, 170-MeV, one, that extrapolate to a value g v l = 0.692(7) 
in the chiral limit, in good agreements with values obtained independently in the meson sector, 
Zy = 0.673(8) and Za = 0.6878(3) [9]. This demonstrates the good chiral and flavor symmetries 
of these ensembles at this relatively low momentum cut off of a = 1.371(8) GeV. 

The isovector axial charge, gA, is also calculated by the corresponding local current, this time 
from averaging its space components. The results, shown in Fig. 2 in the right pane, are noisier 
than the vector charge, especially for the light ensemble which needs more statistics. Fortunately 
the signal for the heavy ensemble is sufficiently good to draw important conclusions as we discuss 
in the following. 
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Figure 2: Left: local-current isovector vector charge gives very clean and accurate signals. We obtain 
estimates of gv = 1.450(4) or 1.447(9), corresponding to Zy = 0.692(7) in the chiral limit, in good agree- 
ment with independent estimates of Zy = 0.673(8) and Za = 0.6878(3) obtained in the meson sector. This 
demonstrates good chiral and flavor symmetries, up to small 0(a 2 ) correction, of the present ensembles even 
at relatively low cutoff of 1.371 GeV. Right: local-current isovector axial charge. 
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Figure 3: Ratio, gA/gv, of the isovector local-current axial and vector charges. Good chiral and flavor 
symmetries of DWF formulation make this ratio naturally renormalized and provide clean and accurate 
signals. We obtain estimates of 1.15(5) and 1.12(10). 



The ratio, gA/gv, of isovector axial and vector charges, shown in Fig. 3, is much less noisy in 
comparison with the axial charge itself. We obtain estimates of gA/gv = 1-15(5) for the heavy, 
250-MeV, and 1.12(10) for the light, 170-MeV, ensembles, respectively. As this ratio is naturally 
renormalized under the good chiral and flavor symmetries of the DWF formalism, up to small 0{a 2 ) 
correction, these values can be compared directly with the experiment of 1.2701(25) [20]: they are 
about 10-% below the experiment. The deficit is almost three-standard-deviation significant for 
the former, heavy ensemble, though the value for the light ensemble needs more statistics to be 
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Figure 4: Non-monotonic dependence of ratio, gA/gv, on m^: it once decreases with nr\ in 1.7-GeV en- 
sembles (red) but returns to an indistinguishable value at a lighter mass in 1.4 GeV ensemble (blue). 

significant. 

Excited-state contamination has been proposed as a possible explanation for this discrepancy 
[10, 11]. The excited-state spectrum depends rather strongly on pion mass: As the separation 
between the ground state and the first excited state is expected to grow with decreasing pion mass, 
such contamination would decrease. Note also the way we calculate the vector charge, gy, protect 
it from the excited-state contamination: conserved charge cannot excite the ground state. In Fig. 
4 we plot our values of the ratio, gA/gv, against the pion mass squared. It is not monotonic: If 
we compare the two points from our previous calculations with Iwasaki action, presented in red in 
the figure, the deficit seems to grow with decreasing pion mass. But the heavier of the two present 
ensembles, presented in blue, brings the value back to less deficit, to a value the difference of 
which is statistically insignificant. This non-monotonic behavior is hard to explain if the significant 
excited-state contamination were present in either of these ensembles, as such contaminations have 
to agree across the ensembles generated with different gauge actions and physics parameters. 

The picture becomes much simpler when we plot the ratio, gA/gv, against the finite-size scal- 
ing parameter, m n L, as is shown in Fig. 3: All the results from our earlier calculations with Iwasaki 
gauge action at cut off momentum of about 1.7 GeV and the present ones seem to align well mono- 
tonically in this plot, with the deficit growing with decreasing m n L. And more importantly the two 
points that share about the same value of m n L ~ 5.8, one from the earlier Iwasaki calculation with 
m n = 420 MeV and L = 2.8 fm, and another from the present calculation with m K = 250 MeV 
and L = 4.6 fm, agrees within each other's statistical errors, 1.19(4) from the former and 1.15(5) 
from the latter, despite their quite different gauge action, cut off momenta, DWF residual mass, 
source-sink separations, quark mass, and spatial volumes. They are very unlikely to agree if one 
or the other or both suffered significant excited-state contamination, as such contamination would 
have to agree with high accuracy despite their quite different gauge actions, cut off momenta, DWF 
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Figure 5: Ratio, gA/gv, of the isovector axial and vector charges, plotted against the dimensionless finite- 
size scaling parameter, m K L. THe dependence is monotonic, and two results at similar values of m K L ~ 5.8, 
1.19(4) from 124 and 1.15(5) from ID, agree with each other, despite very much different m K and lattice 
cuts-off, that significantly alter mass spectrum. 

residual mass, source-sink separations, quark mass, and spatial volumes. 
4. Conclusion 

RBC and UKQCD collaboration continue to calculate nucleon-structure observables using the 
2+ If dynamical DWF ensembles they jointly generate. We found the following from the ongoing 
calculations with ensembles with latticee cut off momentum of a -1 = 1.371(8) GeV, the spatial 
extent L = 4.6 fm and pion mass m % = 250 and 170 MeV: 

1. Nucleon mass estimate is ~ 1.054(7) and 0.984(8) GeV respectively for the heavy and 
light ensembles. 

2. Isovector vector charge is well-behaved, demonstrating good chiral and flavor symmetries 
with DWF at this relatively low cut off momentum. 

3. Isovector axial charge is noisier, yet the calculation is solid. 

4. There is about 10% deficit in gA/gv, with almost 3-standard-deviation significance. 

By comparing the result from the heavy ensemble with our earlier calculation with pion mass 
of about 420 MeV and lattice extent of about 2.8 fm that shares similar value of the finite-size 
scaling parameter, m n L ~ 5.8, we find excellent agreement in gA/gv- 1-15(5) from the present and 
1. 19(4) from the earlier, despite their coming from different gauge actions, lattice cut off momenta, 
source-sink separations, pion mass and lattice volume. Such an agreement is unlikely if there were 
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significant excited-state contamination present in either of the calculations. Thus we conclude 
neither suffers from such a contamination above our statistical error. 

This confirms the scaling in finite-size parameter, m % L. The result can be interpreted as the 
first concrete evidence for the pion cloud surrounding nucleon. With light pion cloud surrounding 
it, nucleon is hardly point-like: at m % = 250 MeV the 4.6 fm linear extent of the present lattice 
is insufficient to contain the axial current carried by the cloud. It grows further toward the lighter 
physical mass of m n ~ 140 MeV. To sufficiently contain the axial charge, spatial lattice extent 
much larger than the 4.6 fm is required. This seems hard to reconcile with the conventional nuclear 
models with point-like nucleons interacting through non-relativistic potential. 

In the talk the status of our structure function study was briefly summarized as well: signals for 
the low moments of structure functions such as isovector quark momentum and helicity fractions 
are noisier than those for form factors, but calculations for these are well under way. 

We are increasing our statistics: double at least by adding extra source coordinates at t = 8, 
24, 40, and 56, and possibly more if necessary, for the heavy ensemble. We plan at least the 
corresponding size of statistics for the light ensemble, using the new AMA technique [21]. We 
seek calculations at physical pion mass, and with appropriate isospin breaking soon afterward. 
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